RU-92-10 
LPTENS-92/11 



(N 



> 

G^ 
in 
o 
in 



Loop Gas Model for Open Strings 



V. A. Kazakov 



! Department of Physics an Astronomy 

Rutgers University, Piscataway, NJ 08855-0849, USA 

and 

Laboratoire de Physique Theorique 
Departement de Physique de I'Ecole Normale Superieure 



^ ! 24 rue Lhomond, 75231 Paris Cedex 05, France 

CN 

I. K. Rostov * 

^ ■ Service de Physique Theorique ^ de Saclay 

CE-Saclay, F-91191 Gif-sur-Yvette, France 

X 

• The open string with one-dimensional target space is formulated in terms of an SOS, 

or loop gas, model on a random surface. We solve an integral equation for the loop 
amplitude with Dirichlet and Neumann boundary conditions imposed on different pieces 
of its boundary. The result is used to calculate the mean values of order and disorder 
operators, to construct the string propagator and find its spectrum of excitations. The 
latter is not sensible neither to the string tension A nor to the mass fj, of the "quarks" 
at the ends of the string. As in the case of closed strings, the SOS formulation allows to 
construct a Feynman diagram technique for the string interaction amplitudes. 

Submitted for publication to: Nuclear Physics B 

SPhT/92-049 

4/91 



* on leave from the Institute for Nuclear Research and Nuclear Energy, 72 Boulevard Tsari- 

gradsko Shosse, 1784 Sofia, Bulgaria 

^ Laboratoire de la Direction des Sciences de la Matiere du Commissariat a I'Energie Atomique 



1. Introduction 

The open string tlieories always attracted a considerable attention of the physicists, 
not only from the point of view of critical strings but also as a possible source of field 
theoretical applications. For example, the idea to formulate the multicolor QCD as a 
theory of noninteracting strings (random surfaces) has been fascinating the minds of some 
theoretical physicists in the 80's. This string theory should involve both closed strings 
describing glueballs and open strings describing the qq bound states (mesons). 

Now, after more than ten years of study, we know how to formulate and solve the 
simplest theory of random surfaces - the (noncritical) closed bosonic string. In order to go 
further, an obligatory exercise to do is to extend the solution to the case of open strings. 

It is clear that the physics of open strings should be more complicated than that of 
closed strings, since it depends on the choice of the boundary conditions at the ends of the 
string. The string amplitudes will depend now on two dimensionful parameters: the string 
tension A coupled to the area of the world sheet (the "bulk cosmological constant" in the 
language of 2d gravity) and the mass fi at the ends of the string coupled to the length of 
the boundary of the world surface (the "boundary cosmological constant"). 

The open bosonic string is well defined for embedding spaces with effective dimension 
("the central charge of the matter fields") — cx) < C < 1; otherwise the vacuum would 
be unstable due to the tachyonic excitation. The field theory of open strings with no 
embedding (C = 0) has been formulated as a random matrix model in [0. This model 
was then solved in the double scaling limit in 0] (see also 0). Further, the theory of 
C = 1 open strings (embedding space IR) was considered as a solution of matrix quantum 
mechanics in [ffl, and 0. Open strings with C = —2 and C = 1 have been also 
considered in [0. 

On the other hand, the noncritical open strings have been studied by means of the 
Liouville theory @, 0, fl^]? IJllIl- '^^^ continuum approach is based on a free field theory 



(the Liouville potential is treated as a perturbation) and therefore cannot be used to 
evaluate the full string interaction amplitudes. This approach is sufficient to study the 
so-called bulk amplitudes which obey the conservation of the Liouville energy. 

In this paper we propose a systematic approach to the open noncritical strings with 
— oo < C < 1 which can be used to find the exact string interaction amplitudes. A very 
convenient framework for this purpose is provided by the loop gas (or SOS-, or heights-) 
model on a random surface [|13| [|14| [jl5| Jl^. The target space in this model is the 



infinite discretized line 2. It is sometimes called loop gas model because the domain walls 
between the regions of constant height form a configuration of nonintersecting loops on 
the world sheet. In the case of Dirichlet boundary condition the points along a connected 
boundary have the same height and therefore the domain walls cannot end at the boundary. 

Below we are going to adapt this model for the case of Neumann boundary conditions 
corresponding to free endpoints of the open string. In this case the domain walls are always 
orthogonal to the edge of the world sheet. This means that the loops are repulsed from 
the boundary but domain walls can approach it at right angle. We will solve the loop 
equation for the amplitude of a disk with a boundary divided into two parts with Dirichlet 
and Neumann boundary conditions respectively. Knowing this amplitude we can further 
calculate the open string propagator and the string interactions following the strategy 



applied in the case of the closed SOS string |]T5|, [ff6[. The eigenstates diagonalizing the 



string propagator are different for different choices of the parameters A and n of the open 
string. However, the diagonalized propagator is universal and is in fact identical to the 
one of the closed SOS string. 

In order to explore the whole range of effective dimensions of the target space —oo< 
C < 1 we will introduce, following the Coulomb gas picture, a distributed background 
momentum ("electric charge") proportional to the curvature of the world sheet metric. 
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The momentum conservation (the electric charge neutrahty) is assured by introducing 
pointhke electric charges at the critical points of the embedding of the world sheet. On 
the lattice this construction has been elaborated in and [1^ 



In the present paper we avoid introducing lattice discretization of the world sheet in 
order to keep closer connection with the continuum theory. A derivation of the basic loop 
equation for a discretized surface is presented in the Appendix. 



2. Definition of the model 

2.1. Coulomb gas picture 

The dynamical fields in the Polyakov formulation of the string path integral [|T9 are the 



position field x{^) and the intrinsic metric GabiO^ a, 6 = 1, 2 of the world sheet (which we 
will denote by We will assume that has the topology of a disk. The boundary dAi 
is divided into 2n pieces as is shown in fig. |l] on which we impose alternatively boundary 
conditions of Dirichlet 

d\\x{^) = t''{C)dax{C) = 0, ta(0 = unit tangent vector (2.1) 

and Neumann 

d±x{^) = n"(^)9ax(^) = 0, Uai^) — unit normal vector (2.2) 

We will denote the Dirichlet boundary by dM.^^^ and the Neumann boundary by dA4^^\ 
Each kind of boundary consists of n connected pieces 

dM^^^ = dM[^^ + ... + dM[i^\ 

c>A^(^) = dMi""^ + ... + dM\^^; (2.3) 



The Dirichlet boundary condition ( |2.1|) is appropriate for the initial anf final string 
states; it describes a boundary which occupies a single point of the embedding space. The 
Neumann boundary condition means that the fiow of energy across the boundary is zero; 
it should be imposed along the edges of the world sheet representing the endpoints of the 
open string. 

The world sheet with n pairs of boundaries describes the interaction of n open strings. 
The corresponding amplitude will depend on the intrinsic geometry of the world sheet only 
through the gauge invariant quantities: the total area of the world sheet 

A= f dA{£)- dA{^) = rf2^v^detG'(0 (2.4) 

and the total lengths of the Dirichlet and Neumann boundaries 

4=/ ^ m). 4=/ m). fc = i,...,n; dm = tamc (2.5) 
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An effective dimension C < 1 can be achieved by introducing a coupling pq ("dis- 
tributed electric charge" in the Coulomb gas terminology) between the field x and the 
intrinsic geometry of the world sheet. The world sheet action then reads 

A[x,Gab] =A'[x,Gab] +A"[x,Gab] 

A'[x, Gab] =£ dAiOG'^'daxiOdbxiO (2.6) 

A"[x, Gab] =tpo \^ [ dAiOxiom + ^ [ diioxiom) 

l-47r Jj^ 2n Jgj^ J 

where -R(^) is the intrinsic Gaussian curvatire at the point ^ E and K{C) is the geodesic 
curvature at the point ^ G dAi. The factor g known as the Coulomb gas coupling constant 
can be elliminated by rescaling x. We fix the normalization of x to have 

9 = l+Po- (2.7) 
The two curvatures are normalized so that the Gauss-Bonnet formula reads 



[ d^^Vdet GR{^) + 2 [ di(^)K{^) = Att (2.^ 

Jm JdM 



'M JdM 

The boundary term in (|2.6|) is introduced in order to be able to satisfy the momentum 



conservation (the "electric charge neutrality" ) . It is clear from the Gauss-Bonnet formula 
( P78| ) that the zero mode x(^) = xq of the a;-field produces only a factor exp(— zpoa^o) and 
can be neutralized by introducing a background momentum — po at some point of the 
boundary. 

Eq. ( pr^ ) defines the standard Coulomb gas description of the C < 1 strings. In this 
paper we propose a modified version of the Coulomb gas approach in which the electric 
charge neutrality is required in a stronger sense. We introduce a system of pointlike electric 
charges associated with the points where the string picture changes. These are the critical 
points of the map — > IR 

rfx(0 = daxdC = (2.9) 

shown in fig. ^. We distinguish four kinds of critical points ^* which will be characterized 
by a weight xiC*) taking values 1, —1, 1/2, —1/2. 

For the critical points in the interiour of the world sheet (cases a , 6 ) we define 

xiC) = sgn det||a„a,x(Olk=e, C eM (2.10) 

For the critical points along the edge of the world sheet (cases c , d ) we define 

Xin = Isgn ajx(e)5=e (2-11) 

The sum over the weights of all critical points gives the Euler characteristics of the world 
sheet 

J2^(e) = X (2.12) 
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Therefore, if we associate with each critical point ^* a charge —poxi^*)^ the electric charge 
neutrality will be fulfilled. The factor 

J^e-^Pox(r) (2.13) 
can be taken into account by adding to the standard action ( p.6|) a second linear term 

A'" = -po [ d\ xiOpiO: PiO = (2.14) 

JM ^* 

The density p(^) can be expressed through the vector field with unit norm n(^) 

P(0 = ^ f d^ix{i)e''h^''dah,d,nd (2.15) 



27r 



M 



MO = ^""fl ... (2.16) 

Consider the functional integral 

Z(A;£„£„a:„i = l,2,...,n) = j [dx{i)][dGabme~^^"''^^'^ (2.17) 

A[x, Gab] = M + A" + A'" (2.18) 

where the integral over intrinsic geometries is restricted to surfaces with fixed area A and 

lengths and of the Dirichlet and Neumann boundaries, correspondingly. The integral 
depends also on the positions Xi of the Dirichlet boundaries in the embedding space. The 
interaction amplitude of n open strings with momenta and lengths is 

defined by the Laplace transform 

v{pkJk;k = l,...,n) = 

poo n „oo _ _ _ (2 19) 

/ dAe-^^^W rf4e-^''=+*P'="'=Z(A;£„4,x„i = l,2,...,n) 
Jo ^^1^0 

Here x^ denotes the position of the /c-th Dirichlet boundary dM.\^\ The string tension A 
coupled to the area of the world sheet is called sometimes cosmo logical constant, since this 
functional integral can be also considered as the partition function for two-dimensional 
quantum gravity. Similarly, the mass p of the ends of the string can be also interpreted 
as a boundary cosmological constant since it is coupled to the length of the Neumann 
boundary. 



^ Strictly speaking, this term is not linear in x because of the charge density p(^) depending on 
the embedding ^ Note also that each connected Dirichlet boundary contributes a factor 

— Po if it is a closed loop and — po/2 if it is an open interval. 
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The amplitude ( p.l9| ) is nonzero only if the sum of all momenta is equal to the back- 
ground momentum (1 — n/2)pQ 

Tl 

+ = (1- -)po (2.20) 

It is convenient to introduce the variables Zk dual to the lengths ik and consider the Laplace 
image of ( |2.19| ) 

v{pi,Zi;...;pr,,Zr,)= dh... e-^i^i--^-^"^(pi,£i;...;p,,C) (2.21) 

Jo Jo 

The presence of the background momentum po diminishes the effective dimension of 
the embedding space (the conformal anomaly due to the matter field) from one to 

C = l-6pl/g = l-6{g-l)yg (2.22) 
and restricts the spectrum of allowed momenta to 

p = |po, keT (2.23) 

The local operators in the theory are those creating microscopic closed and open 
strings. Sometimes they are called bulk and boundary operators [|lO|. The spectrum of 
the bulk operators can be fixed with self-consistency arguments known as David-Distler- 
Kawai analysis [|19| , based on the assumption that at distances large compared to the cutoff 
but small compared to the size of the world sheet, the fiuctuations of the metric Gab are 
described by a gaussian field. 

Below we present a sketch of these arguments mainly to help the reader to become 
familiar with our normalization which is not the standard one used in the string theory. 

After introducing a conformal gauge 

Ga6 = e2-'^(«)G°,(0 (2.24) 

where G^j^ is some fiducial metric, and taking account of the conformal anomaly we arrive 
at an effective action depending on a two-component gaussian field (x, (p): 

^Liouville[x, <P] = ^ [ d^JdetG^iOig G'^'^'idaXiOdtxiO , ^ 

J ^ (2.25) 

- darndbHO) + i^po xio - eom)R\o] 

The vertex operator creating a momentum p, dressed by the fiuctuations of the metric 

V(p^,)(0 = e*(p-P")^(«)e-(^(p)-^o)^(«) (2.26) 

In particular, the puncture operator V = —d/dA which marks a point on a surface is 
represented by 

p(^) = Q-i'^iPo)-^o)'P{0 — g2i/<^(?) (2.27) 



Here we write explicitely the compensating charge —po associated with the puncture 
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The condition of the absence of conformal anomaly yields 

Ctot = C(,) + C(^) - 26 = 1 - 6^ + 1 + 6^ - 26 = ^el-pl = Ag (2.28) 

We choose the positive solution thus fixing a positive direction in the space 

eo = (7 + l, Po = g-l (2.29) 
The condition that the conformal dimension of the operator ( pj.26| ) is one 



= 



(2.30) 



combined with (|2.28|) leads to the relation 

e{pf-p^ = Q (2.31) 

which can be interpreted as a mass-shell condition for the 2-momentum (p, e) . All physical 
operators correspond to positive Liouville energies 

e{p) = \p\ (2.32) 

The Liouville charge of the identity operator equals to eo — e(po) = 2z/ where 

v = ]^{g + l-\g-l\) (2.33) 

The gravitational dimensions of the vertex operators coherent with the background mo- 
mentum pq — \g — 1\ are 



5., = 1 - li^IM^ = I' -^-l-l^- I (2.34) 



eo - e(Prs) _ 


r — gs 




\9-^ 


eo - e(po) 


9 + 1 




g-l\ 



Finally, the string susceptibility exponent 7str giving the dimension of the string interaction 
constant is equal to 

e(Po) 2|£^-1| 
lstr = = , , , I 77; z/(2 - 7str-) = eo (2.35) 

The above arguments can be easily generalized to the boundary operators How- 
ever, as it has been noticed in |]T^-[|T6[, the semiclassical analysis is not always applicable 
at the boundary. 

2.2. Formulation as an SOS model on the world sheet 

Let us now try to find a link between this continuous formulation of the path integral 
and the so called SOS model in which the x field is restricted to take only discrete values 
(heights) x/tt G 2. At large distances the configurations of such field should look as 
continuous; this is achieved by the condition that the x field can jump only with a step 
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±7r. The domain walls separating the domains on the world surface where x takes constant 
value form a pattern of nonintersecting lines. In the case of a surface without boundary 
all these lines should be closed loops. Across each domain wall the heigth x jumps by n : 
X X ± n. The sign can be taken into account by assigning an orientation to the domain 
wall. 

If we consider a Dirichlet boundary, the above picture holds unchanged. Since the 
height X is not changing along the boundary, the whole boundary belongs to a single 
domain. Note however that the loops are allowed to touch the boundary and this should 



be taken into account when writing the loop equations for the closed string [p!2|. 

On the other hand, in presence of a Neumann boundary the above geometric picture 
changes drastically. The domain walls are not only loops but also lines ending at the 
boundary (fig. The condition that the normal derivative of the x field is zero in the 
vicinity of the boundary means that all these lines meet the boundary at right angle. In 
addition, the closed loops are not allowed to approach the boundary. 

The integration over the x field can be replaced by a sum over all loop configurations 
on the world sheet and a subsequent sum over all allowed values of x in the domains 
bounded by these loops: 

l[dx{0...-^ J2 E - (2.36) 

^ loop configurations x(domains) 

Suppose that the integral over the world-sheet intrinsic geometries is regularized, say, 
by a discretization using planar graphs. Then ( |2.36| ) can serve as a microscopic definition 
of the string path integral. 

Let us "derive" the Boltzmann weights of the domain-wall configurations from the 
continuum action ( |2.18| ). We can imagine that the SOS configuration is regularized so 
that the map Ai — > 7r2 is obtained as a limit of a smooth map — > IR. Then the 
contribution of the last term in ( |2.18| ) comes only from the vicinity of the domain walls. 
Notice that d<f>{^) = n(^) x dan{^)d^°' is the infinitesimal angle swept by the unit vector 
UaiO along the interval d^. 

Let us consider an SOS configuration of the field x(^) described by a system of domains 
and domain walls and evaluate the action ( p.l8| ). 

The term daxd°'x in the integrand is just a square of the invariant gradient of the x- 
field, which is zero everywhere except of the domain walls, where it is an (infinite) positive 
constant along the wall. Being integrated over the world surface it yields the total length 
of the domain walls times a (cut-off dependent) positive factor. Thus its contribution to 
the action is 

A' = Ko [ dm = ^o4otai (2.37) 

J domain walls 

where dl{^) is the length element along a domain wall and ^totai is the total length of the 

domain walls on the world sheet. 

tNow let us demonstrate that he contribution of the last two terms in ( p.6| ) depends 

only on the topology of the configuration of domains and domain walls. Let us consider 
a world sheet with the topology of a disk and a system of domain walls separating the 
domains Vi,V2i ... on it. The domain is bounded by domain walls (loops and 

open lines) and pieces of the boundary of the world sheet. Two neighbour domains can 
be separated either by a closed loop or by an open line. For each two domains having 
a common boundary one of them is surrounded by the other; therefore the system of 
domains has a tree-like structure. The boundary of each domain consists of connected 
components. The Cjt — 1 internal boundaries are all closed loops. The externel boundary is 
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made out of open lines separated by pieces of the boundary of the world sheet, = 

means that the external boundary is also a closed loop. 

Now we can express the contribution of the term in ( |2.6| ) proportional to the gaussian 
curvature in terms of the heights Xk of the domains Vk and the numbers Ck^Uk character- 
izing the topology of the domain wall configuration. Applying the Gauss-Bonnet formula 
to each domain domain Vk we find 



A"-z^ f diiOxiOKiC) 

= [ dAiOxiom 



M 

(2.38) 



V X J2 ( / dl{i)k{^) + TTTifc) + 47r(/ifc + Cfc - 2 



k 

■ Po 



^47r 

k 

where K{^)di{^) is the infinitesimal angle swept by the normal vector na(0 along the 
boundary dVk, and hk denotes the enclosed genus handles) in the domain Vk- Since 
we are considering the topology of the disk, hk = 0. The boundary integral is understood 
as a sum of the integrals over the smooth pieces of the boundary. The last term on the 
r.h.s. is due to the most external connected component of the boundary dVk of the domain 
Vk having 2nk edges with angle 7r/2; their contribution to the global geodesic curvature is 

TTHk. 

Now let us consider the third term A'" . After integrating by parts the integrand 
in ( p.l5|) turns to dx{$,) /\d(p{$,). It is easy to see that the contribution of each domain 



wall is ^|^(a;I■ight — a:ieft)v'giobai where </?giobai is the angle swept by the normal vector n(^) 
along the domain wall (it is equal to the integral of the geodesic curvature). Adding the 
contributions of all domain walls we find 



JdM 47r jj^ 



A'" + 1^1 dmxiom) = ^^l rf"^a„x(^)£«Wa;,nrf 

(2.39) 



.Po 



Collecting the three terms ( |2.37D , ( |2.38| ) and ( |2.39| ) we arrive at the following action for 
given loop configuration 

A = A' + A" + A'" = Koitot + ipo J2^k[2-Ck- \nk\ (2.40) 

k 

Thus the Boltzmann weight e-^ of each domain wall configuration depends only on its 
topology and the total length of the loops. The action ( |2.40D can be simplified further by 
performing the sum over the heigths Xk of the domains T>k- We have to calculate the sum 

Q ^ ^gipoEfc^'=(2-'^fc-5"0 (2.41) 
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The calculation is performed in the same way as in the case of the closed string [12| |16 . 
We will exploit the fact that the system of domain loops and open lines on the disk has 
a tree-like structure. Let us start with a domain T>k on the top of the tree, i. e., a 
simply connected one. It is represented by a vertex with a single line (tadpole) of the 
corresponding graph. Consider first the case when nfc = 0, Cfc = 1 when the boundary is a 
closed loop. Then the sum over xinside — Xk yields 



g^Poxin.idc ^ 2 cos(7rpo)e'P°^°"*=''^« (2.42) 

■^inside -^outside 

But Xoutside IS the X Coordinate of the surrounding domain. Therefore the result of the 
summation is a factor 2 cos(7rpo) and a reduction by one of the number of connected 
boundaries ( c — > c — 1) of the surrounding domain. Proceeding in the same way we can 
eliminate one by one all loops until we arrive at a configuration (a "rainbow diagram") 
containing only open lines ending at the boundary of the world surface. Each domain Vk 
is characterized by the number rik of the domain walls along its boundary [ck = 1). This 
configuration has again a structure of a tree and we can sum over x as before starting with 
the domains on the top of the tree, i.e., these whose boundary is formed by a single line 
(cfc = 1, n/c = 1). The sum over the x coordinate of such domain yields a factor 2 cos(^7rpo) 
and eliminates the term associated with its boundary 



g^ipoxinsidc = 2cos(-7rpo)e^'^"''°"*='''" (2.43) 
_ _i_ 2 

^inside — ^outside itTT 



E 



After repeating this procedure several times we eliminate all domain lines. Thus the sum 
over the embeddings produces the following weight of each configuration of domain walls 

/ \^ loops / ][ \ #open lines 

0= (^2cos(7rpo)j (^2cos(-7rpo)j (2.44) 

The sum over the last coordinate yields an infinite factor which is the volume of the 
embedding space. 

Summarizing, we arrived at a modified loop gas model on the random surface. Its 
partition function is a sum over configurations of nonintersecting loops and open lines 
ending at the boundary 

^= E E 

surfaces loop configurations ^2 
/ \ #open lines / \ #loops 

g-2KoAot 1 2 cos(7rpo/2) j f 2 cos(7rpo) j 

The construction of the generalized loop gas can be made more explicit by discretiz- 
ing the measure over random surfaces as prescribed in |]T^. The only difference is that 



the curvature is concentrated at the sites of the lattice and the Gauss-Bonnet theorem 
degenerates to the Euler formula. 



3. Loop equations for the open string 

In order to exploit the definition (|2.45| ) we have to give a meaning of the functional 
integral over surfaces. It is convenient to take the two sums in ( |2.45| ) in the opposite order: 
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first to fix the topology of the configuration and the lengths of all lines, and then perform 
the sum over the geometries of the connected pieces of the surface (the "windows" ) . Each 
"window" contributes a factor depending only on the length of its boundary. Finally we 
integrate over the lengths of the lines and sum over all topologies. This sum can be most 
easily performed using equations of Dyson-Schwinger type [0. Below we will use the 
continuum formulation of the Dyson-Schwinger equations proposed in [pO |. 



In order to obtain a closed loop equation we have to consider a disk with only one 

pair of Dirichlet and Neumann boundaries with lengths I and I. It seems that the only 
consistent way to avoid loops touching the Neumann boundary is to have an open line end 
at each point. If we are using a lattice regularization, this means that a line is ending at 
the middle of each bond forming the Neumann boundary (see Appendix A). 

Let V{i., i) be the partition function of the disk with such mixed boundary condition^. 



It is related to the functional integral (|2.17|) with n = 1 b 



ViiJ)= (iAe-^^Z(yl; £, £, x) (3.1) 
Jo 

An infinitesimal deformation of the Neuman boundary at its endpoint (one of the points 
separating the two boundary conditions) singles out the line starting from this point which 
splits the world surface into two pieces. The loop equation follows from the geometrical 
decomposition of the disk shown in fig. § 

(9 _ ~ 

—^V{i,i) =2 cos{7ipo/2) di' di'e-^^°^'v{i',i') V{i-i',i + i') (3.2) 
di Jo Jo 

Eq. ( p.2| ) has a clear geometrical meaning. It sums up the rainbow diagrams with an 
additional structure: the space between its lines is occupied by surfaces with loops. Note 
that this loop equation determines only the dependence on L; therefore it has to be com- 
plemented with another equation specifying the dynamics of closed loops. The missing 
information can be supplied by fixing W{i) = V{0,i) which is exactly the partition func- 



tion of a disk with Dirichlet boundary conditions. It satisfies a loop equation |]T2| - |T^ of 
the type 



-U'[^^W{i)= j di'W{i')W{i-i') 



+ 2cos(7rpo) / di'e-^^^^'W{i')W{i + i') 
Jo 



(3.3) 



where U' is some local (differential) operator describing an infinitesimal deformation 

of the boundary of the disk. 

As usual, in order to turn the convolution in ( |3.2| ) into a product, we introduce the 
Laplace transform 

V{T,t)= dL die-^^-^^V{i,i) (3.4) 
Jo Jo 



^ We denote this amplitude by V saving the letter v for the corresponding renormalized 
amplitude 

^ When n = 1 the open string amplitude does not depend on the position x of its only Dirichlet 
boundary 
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and eq. (|3.2|) turns to 



TV{T, t) - W{t) = ^'''^^Pom j ^v{T, t')V{T, 2Ko - t') (3.5) 
where T plays the role of bare mass of the "quarks" at the ends of the open string and 

/"OO /"OO 

W{t) = / dTV{T, t)= di e-^^W{i) (3.6) 
Jo Jo 

The contour integral in ( |3.4| ) goes around the singularities of W{T,t) and leaves outsides 
the singularities of W{T, 2Kq — t). Similarly, the Laplace transform of (|3.3|) reads 

W{t)^ = l^.f j^W{t'f[U'{t') - 2 cos(7rpo)Vr(2i^o - t')] (3.7) 

All these loop equations can be derived in a rigorous way starting from the lattice version 
of the model (see Appendix A). 

The loop amplitude V{1, T) can be considered as the classical background field in an 
open string field theory. It satisfies a mean-field type equation which is equivalent to eq. 
( p.2|) . This equation is derived by cutting the world sheet along the most internal open 
lines as shown in fig. |5|. In this way the amplitude V can be^^xpressed as an integral of 
the product of a VF-amplitude and a number of V amplitudesQ 



V{t,T) 

E / die-'' n ^e-2^o^^2cos(7rpo/2)y(4,T) W{l + h + 

n=0-^0 k=l \ J 

dt' Im#(2i^o - 1') 



/ 

J M 



7r{t + t') 1- 2cos{npo/2)V (T,t')/T 



It is known 



14 16 21 



that depending on the explicit form of the operator U{d/di) 
one can achieve different critical regimes at the critical temperature K* of the loop gas on 
the random surface. Here we will consider in details the so-called dense phase corresponding 
to the simplest choice U'{d/dt) = djdl. In this phase the loops fill the world surface 
densely, without leaving space between them. One of the peculiarities of the dense phase is 
that the fractal dimension of the Dirichlet boundary is larger than one: Xjv = 1/(1 — |po I) = 
1/(7. The dilute phase of the loop gas corresponds to (multi) critical potential U [1T6| [pT| . 
The potential is tuned so that the area of the world surface not occupied by loops also 
diverges. The equation for the loop amplitudes is the same for both phases but the scaling 
of the cosmo logical constant is different. In the dilute phase the fractal dimension of the 
Dirichlet boundary is 1/z^ = 1. In the Coulomb gas picture the dense and dilute phases 
are related by a duality transformation g ^/g- 



^ By V{l,T) we denote the Laplace image of (3J_) w.r. to £ ; it depends on T through the 
factor exp (# open lines) 
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4. Solution of the loop equation in the scaling limit 



We will follow the method worked out in [|T2[ for solving eq. (|3.3| ) directly in the 
continuum limit, and apply it to our master equation (|3.5|). 

Let us first recall the solution of eq. (|3.3| ). We expect that W{t) has a cut t^, tn along 
the real axis of the t-plane (on the first sheet of its Riemann surface). The contour of 
integration in ( |3.7| ) goes around this cut. If we replace in ( |3.7D t with 2Ko — t the integrand 
will not change, but the contour of integration will envelop the cut [2Kq — t^^, 2Kq — t^] of 
the function W{2Ko — t). Therefore, adding these two equations, we integrate along both 
contours which form together a contour surrounding all singularities of the integrand. 
Applying the Cauchy theorem we find the following functional equation for W{t) (we 
consider the simplest differential operator U'{t) = t) 

W{tf + W{2Kq -tf + 2 cos{TipQ)W{t)W{2Ko - t) 

^tW{t) + {2KQ-t)W{2KQ-t)-2 (4.1) 



Taking the imaginary part of (|4.1| ) and knowing that lmW{t) ^ along the cut, we arrive 
at a linear Cauchy- Riemann problem: 

ReW{t) + cos(7rpo)Inil¥(2Ko -t)= t/2, t G [^l, Ir] 

ImW^(t)=0, t^[tL.tR] (4.2) 

If we take eq. ( |4.1| ) at the symmetry point t = Kq we obtain 

W(Ko) = ^ (4.3) 

Ko + V^o'-4(l + cos(7rpo)) 

In the dense phase the temperature 2Ko of the loop gas is also the bare cosmological 
constant since the total length of the loops is equal to the area of the surface. The 
singularity of eq. ( [4.3|) gives its critical value 



Kq = 2^1 + cosTrpo = 2^2 sin(7r£//2) (4.4) 

At that point the two cuts touch each other: 

t^R = 2K* -t*R = K* (4.5) 

In order to explore the vicinity of the critical point we blow up, as usual, the infinites- 
imal vicinity of the point t = Kq by introducing a cutoff parameter a playing the role of 
elementary length along the boundary 

t = K^ + az, tR = K^- aM (4.6) 

The parameter z is coupled to the renormalized length of the boundary (a boundary 
cosmological constant) and M is the renormalized position of the cut. Note that the 
characteristic length of a loop grows near the critical point as {Ma)~^. 
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Since the singular part of the loop amplitude behaves for M = as 

W{t)r^{azr, g = l-\po\ (4.7) 

we define the scaling part of W as 

W{t) = W* + Aga<^w{z) (4.8) 

where W* is the critical value of W at t = Kq , Kq = Kq and Ag is a constant factor 
depending on the normalization of w. 

Finally, we introduce the renormalized cosmological constant 

Ko = K* + Bga^'^A (4.9) 

where Bg is an appropriate constant and u has the meaning of the inverse fractal dimension 
of the Dirichlet boundary, if the dimension of the world sheet is assumed to be 2. Since 
A is the only parameter in the theory, we expect that M^*^ ~ A. To determine u we note 



that from (|3)-(|J) (for l/2<g <1)U 

W(Ko) =W* + Ayw(0) 

. ; , (4.10) 

= VT* - (2/K*)3/2/b;A 

and, sinse w{0) ~ M^, we find v = g. 

The renormalized loop amplitude w{z) has a cut — oo < z < —M and satisfies the 
following equation which is a direct consequence of 



(4.11) 



Rew(2;) — cos['ng)w[—z) =0, z < —M 
Imw(2) =0, z> -M 

If we parametrize z by means of a new variable r 

^ = Mcoshr (4.12) 
the reflection z —>■ —z ± iO corresponds to r —>■ r ± in and ( |4.11J ) is replaced by 

^^ind/dr ^ ^-^n^/^T _ 3 cOs{7r g)]w{z) =0, g = 1 - \po\ (4.13) 

with an evident solution 

cosh((7r) 



w{z) = -M^ 



cos((77r/2) 



(z + Vz^ - M^)9 + {z + Vz^ - M2)s 
2cos(£?7r/2) 



(4.14) 



^ In this interval Kq = Kq up to terms of higher than hnear order in the cutoff a 
^ This normahzation corresponds to Ag = 2^3/2-1/2^^ _ ^)9-i[sin(7r5r/2)]"^"^ 
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We have normalized the solution so to have 

w{0) = -MS (4.15) 

Then by ( P^) , with Bg = [K^/2]^Al , the relation between M and A is just 

A = M^s, 0<g<l (4.16) 



The function ( 4.14|) has a cut [— oo, — M] on its physical sheet whereas the cut [M, oo] 



appears only on the next sheets. 

In the same way we can analize the dilute phase of the loop gas {g > 1). We would 
obtain the same expression ( [4.14| ) for the loop amplitude but the scaling of the cosmological 



constant will be different: A = M [|T6[ [^. The scaling of A in both phases of the loop 
gas is determined by the dimension do = ^/t^ oi the Dirichlet boundary, with u given by 

(EH) 

A = M2-,z.= |?' (4.17) 



Let us now consider the loop equation (|3.5|) for the open string. We choose to work 
in the dense phase, but all calculations can be easily extended to both phases of the loop 
gas. Again, after symmetrization w.r.t. the reflection t — > 2Ko — t 

T[V{T, t) + y (T, 2Kq - t)] = W{t) + Wi2Ko - t) 

. . (4.18) 
+ 2 cos(7rpo/2)y(T, t)V{T, 2Ko - t) 

We have assumed that V{t) has the same cut as W(t). 

In order to determine the critical value of T we consider eq. ( [4.18| ) at the point t = Kq 
where it becomes algebraic 

2TV{T,Ko) = 2#(Ko) + 2cos(7rpo/2)y2(T,Ko) (4.19) 

Using ( |4.1C1| ) and dropping all powers of the cutoff a higher than a^/^ we write its solution 
in the vicinity of the critical point as 



T - a/T2 - 2W{Ko)2 cos(npo/2 
V{T,Ko) = — 



2cos(7rpo/2) (4.20) 



V* - aS/^^/2cos{7vpo/2)JAg{^I + 2^1) 



where 



2W* 



is the critical value of the open string amplitude and 

2aSfi = Ag[T''-T^] (4.22) 
14 



is the renormalized mass at the endpoints of the open string (the parameter coupled to 
the length of the Dirichlet boundary) . 

As before, we retain in the scaling limit only the singular part of W{t): 



V{T,t) = V* + a' 



A, 



2 cos(7rpo/2) 



(4.23) 



Comparing (11:231) and (14201) we find 



2(^ + VA) 



(4.24) 



Then, throwing away the higher order terms we obtain from ( |4.18| ) a = g/2 and 

w{z) + w{—z) + v{fx, z)v{fx, —z) = 2fj, (4-25) 



At 2; = this equation reproduces (|4.24| ). 

This equation is compatible with the integral equation 
deed, introducing the scaling variables according to 
integral equation 



(3.8) in the scaling limit. In- 
(|3) and O) we find the 



v{z, fl) 



dzi w{-zi) 
-im- 



M n{z + zi) v{zi,y) 



(4.26) 



Taking the imaginary part of ( [4.26| ) along the cut we find 

Imw(z) 



lm.v{z, n) 



v{z,ij) 



z < -M 



(4.27) 



which gives the imaginary part of (|4.25| ) . 

Exactly the same equation can be obtained for the dilute phase of the loop gas on 
the world sheet which corresponds to the choice g > 1. The only difference is that the 
cosmological constant is replaced by A = M^. All further arguments are valid for both 
regimes. 

If we parametrize z by ( [4.12| ) and /U by 



/i = cosh{ga)] v{fj,, z) = f (cr, r) 



(4.28) 



eq. ( [4.25| ) becomes 



/ . N / N , ( cosh[gf(r + in)] + cosh(gfr) „ , , 

y cos(5f7r/2) 

After shifting r to t + in /2 , eq. (|4.29| ) becomes 

v{t + in/2)v{T - in/2) = AM^ cosh[-(r + a)] cosh[-(r - a) 

2 2 



(4.29) 



(4.30) 
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In the limit A = 0, ;U = the solution of ( [4.3q) is 

v{z) = {2zy/^ (4.31) 
For for nonzero A and ^ but for some particular values of a 

a = ±in/2,±in/2±in/g (4.32) 
eq. ( |4.3C1| ) has a solution in elementary functions 

v(±iTv/2, t) = COSlli OT / Z] 

(4.33) 

v{±in/2 ± in/g, r) = -2M^/^ smh{gr/2) 

In order to solve eq. (|4.28|) in the general case let us take the logarithm of both sides 
to obtain a linear equation on 

u((T, r) = logf (cr, r) (4.34) 

of the form 

(e^f^ +e"'^^)w(a,r) = log[2 cosh(£^(r + a)/2)] + log[2 cosh(£^(r - a)/2)] (4.35) 

It is easy to solve it by performing a Fourier transform which gives an integral representa- 
tion for w(a, r) 

w(a, r) = u{t, a) = f{T + a) + f{T - a) (4.36) 

1 duj e^"^^ 
f{r) = f{-r) = -- J^^ ~cosh{7rLj/2)smh{7ru/g) ^^'^'^^ 

The ambiguity due to the singularity at w = is lifted by imposing the condition /(r) 
g\T\/4 when r ^ oo. By deforming the contour of integration and applying the Cauchy 
theorem we can write the integral ( [1.37 ) as the following formal series which makes sense 
for Rer positive 

r, 1 °° ( ^—gnr °° / — 1 _ — (2A:-1)t 

■'^ ^ 4 2^ n cos{TTgn/2) ^ 2A; - 1 sin[(2A; - l)7r^] ^ ' ^ 

When g is rational, this series for df /dr can be easily summed up. 
Consider for example the case when g = p/q with p, q co-primes and p even. In this 
case q is automatically odd. 

Representing the summation indices in ( [4.38| ) as 



n = 2qN + n, = 0, 1, 2, n = 1, 2, 2q 
k = pK + k, K = 0,1,2,... -,1 = 1,2,..., p 
we arrive at the following expression 

dfir) 1 ( , „ ^n-v/2 cosh[^(g - n)r] 



(4.39) 



g cosh(pr) +9^{ 



dr 4sinh(pr) y cos[g{q — n)^ /2\ 



I 2 Vf 1^ ^osh[(p-H -2k)T] 
^ sin[(p+l-2^)7r/(7] 

fc=i 
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(4.40) 



In the simplest case g = 2 (p = 2,^ = 1) this expression reproduces the result 



fir) = log(cosh ^); ^ = lth(r/2) (4.41) 

which can be easily obtained directly from the integral ([4.37| ). 

If = p/q with p odd, some of the coefficients in this series become infinite. This 
happens when ng = 2k — 1. It is easy to see that the divergent coefficients appear in pairs 
and the contribution of each such pair is finite: 



T 



II r » ^ ^g-(2A:-l)T (4 42) 

n 2cos(7r^n/2) 2k - 1 sin[(2A; - l)7r/^] n n ^' ^ 

The r.h.s. represents the limit of the l.h.s. when g —>■ {2k — l)/n. For example, for g = 1 
we obtain: 

/(r) = - log(cosh r) + -arctge"" - — 1 ^..^ (4-43) 
4 TT 27r [2n + 1) 

This result is already inexpressible in terms of elementary functions, unlike the formula 
for the derivative df /dr. 

Let us note that our disk amplitude f (a, r) being represented in the integral form ( ^.37| ) 
is remarkably similar to the 5'-matrix of thej9(n) -vector model with n = —2 cos{'Kg) on the 
regular lattice presented in the paper [^. la. Our r-parameter corresponds to the rapidity 



parameter in the two-particle S'-matrix. Eq. (4.30) is analogous to the unitarity condition 



on the S'-matrix. This S- matrix was first calculated in |2^] in terms of an infinite product 
of gamma functions, which we can use for our amplitude as well. Expanding cosh ^ in 
( [1.37|) in the exponents and performing the integration we obtain: 

0^ _g2/(r) _ ^ TT U2 +^ 4 +2^zJM2+^ 4 2.^) ^ ^2+^ 4 ) u 44\ 

^ ' J-l 1 , Mi , Z9 )r(i + - ^) r2(i + Q^) ^ ' 

The exponential factor in front of the product depends on how we treat the singularity at 
= in this integral. It is defined through the asymptotics ( |4.31j ) of f (r). 

It is not clear whether this coincidence is accidential or it refiects some deep relation- 
ship between the 0(?T.)-vector field in the fiat and fiuctuating metric of two-dimensional 
space, respectivly. May be the representation of the model in the fiat space in terms of the 



effective Sine-Gordon theory presented in [23| can shed some light on this strange fact. 



5. Boundary operators 

It is very convenient to regard the loop amplitudes as expectation values of operators 
creating boundaries on the world sheet. For this purpose we are going to introduce the 
following notations. Denote by 0{£) the operator creating a closed Dirichlet boundary of 
length £ on the world sheet. The expectation value of this operator is nothing but loop 
amplitude (|4.14| ) 

Q roo 

u,{z) = - — {0{z)), 0{z)=J^ die-^'Oii) (5.1) 



notice a misprint there:the factor 1 /k was missing there in the integral 
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(The derivative comes from the loop amphtude being defined with a marked point on it.) 

Similarly, by 0{i) we denote the operator creating closed Neumann boundary of length 
i. Its expectation value is the loop amplitude with Neuman boundary condition 

w{^) = - {Oif,)), 0{ii)= Me-'^0{1) (5.2) 

Once a Dirichlet boundary exists, one can define a boundary operator C{t) creating 

an open Neumann boundary of length £ at some point. In a similar manner we define the 
operator C{t) creating open Dirichlet boundary of length I at some point on the Neuman 
boundary. By construction 

C{e)0{£) = C{£)d{i) (5.3) 

The disk amplitude v{z, /i) with Dirichlet-Neumann boundary conditions is the expectation 
value of any of the products (| 



v{z,fi) = {C{fi)0{z)) = {C{^,)0{z)) (5.4) 

Here we used the notations 

C{z)= / die-''^C{i), C{fi)= / die-^^Cii) (5.5) 
Jo Jo 

The operator C{i) (resp. C{i) ) creating open Dirichlet (resp. Neumann) boundary 
can be expanded as an infinite series of local boundary operators in the same way as the 
loop operator in the closed string is expanded as a series of operators creating microscopic 
loops [^. (The boundary operators in the framework of the Liouville theory have been 
studied in 0. ) 

Consider first the limit of large /i corresponding to small Neumann boundary and 
expand v{fx, z) in negative powers of ii 



cos((77r/2) sm(7r/(7) 



-k—n/g 
k,n=0 



(5.6) 



~ ~ M9cosh{gT) ~ Mcoshr 

Co,o = 1, Ci^o^ ——=w[z), Co,i = z,... (5.7) 

cos(7r(7/2) sm(7r/(7) 



with 



The coefficients in ( |5.6| ) can be interpreted as expectation values of local boundary opera- 
tors H 

Ck,n{^,z) = {Ck,nO{z)) (5.8) 



^ The operators Co,o,Co,i,.-- are not, strictly speaking, local operators. They are "boundary 
operators" for the Dirichlet boundary 
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so that ( |5.6| ) would imply the following expansion of the operator C(/x) creating Neuman 
boundary as an infinite series of local scaling operators 



C» = 5^4,„A*^-'=-t (5.9) 

k,n 



The leading nontrivial coefficient in ( |5.6|) is just the amplitude 'w{z) of the closed string 
and the corresponding boundary operator is 

Ci,o = (5.10) 

Let us define the dimension 5 of a local boundary operator C acting at a point of the 
Dirichlet boundary. The mean value {CO{z)) = -F(A, z, fx) of such an operator is assumed 
to have the following scaling property 

where 

a = {2 - rstr)x - doil - S) (5.12) 

Note that we measure the dimension of the operator C in units of dimension of the 
Dirichlet boundary and not the world sheet. Otherwise we would have an additional factor 
of dD/2. Let us recall that the dimensions of the Dirichlet boundary is do = 1/u and the 
dimension of the Neumann boundary is = l/u = g/u if the dimension of the world 
sheet is 2. The term d^ (resp. d^ ) comes from the fact that marking a point on the 
boundary breaks the cyclic symmetry and produces a factor of length. In the case of the 

topology of a disk (Euler characteristic (x = 1) ), eq. (|2.35|) yields a — and the 



dimension of the operator C is related to a hj S = a — g/u. 

Now let us examine the mean values Ck^n- It is easy to check that they satisfy 
the scaling ( |5.11| ) with a = {kg + n)/u. Therefore the (boundary) dimensions of the 
corresponding local operators are 

5k,n = {k-l)g + n (5.13) 

The dimension of the operator Ci^ is zero, as expected. 

Let us consider the opposite limit z ^ oo corresponding to small Dirichlet boundary. 
The corresponding expansion of v is obtained from ( [4.36| ) and ( [1.38| ) with a > r > From 
( [1.36|) and ( [1.38|) , assuming that r > cr > 0, we find 



;(r, a) = M^/V^/^fl + ^"^^(^a) ^ 



cos((77r/2) sm{'K/g) 



cosh a _ , 
e-" + ... 



-k)g- 



(5.14) 



/c,n=0 



with 



Mgcosh(ffa) Mcoshg 

Co,o = 1, Ci,o = -, TTTT /U, Co,i = = w[iJ,), ... (5.15) 

cos{TTg/2) sm{n/g) 
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Analogously to the previous case we define the dimension 5 of a local operator at some 
point on the Neumann boundary by the scaling properties of {CO^y)) = -F(A, /i). In this 
case the power a in ( |5.11j ) is related to the dimension of C by 



Q! =(2 - 7st^)x - (i7v(l - 5) 

(5.16) 

={^+gS)/u 

We consider the coefficients Ck,n as mean values of microscopic operators Ok,n at the 
Neumann boundary which implies the expansion 

C{z) = J2Ck,nZ^^'--'^-^ (5.17) 



The coefficient Ck n obeys the scaling ( p.ll| ) with a = {gk + n) /z/ and we find by 

ft — 1 

5k,u = k + (5.18) 

9 

The identity operator is Oi^ = —d/dz. Its expectation value Co,i gives the loop amplitude 
with Neumann boundary condition 



2sm(7r/^) 

Let us make the following remark. The duality transformation of the functional inte- 
gral with the gaussian action (|2.6| ) leads to a similar action but with g replaced hy 1/g and 
the Dirichlet and Neumann boundary conditions exchanged. The symmetry of the func- 
tion u{a, r) is a manifestation of this duality symmetry. In this sense the dense {g < 1) 
and the diluted {g > 1) phases of the loop gas are dual to each other. In the dense phase 
the Neumann boundary has the classical dimension = g/u = 1 while the Dirichlet 
boundary has anomalous dimension (i£) = l/z^ = l/(7>l. In the dilute phase the Dirichlet 
boundary has classical dimension do = l/z^ = 1 while the Neumann boundary has anoma- 
lous dimension d^ = g/i^ = g > ^- The loop amplitude with Neumann boundary condition 
is related to that with Dirichlet boundary condition by z <-> ^, (7 ^ l/^r. Therefore the 
quasiclassical treatment (see, for example [^^) is applicable for the Neumann boundary 
when g < 1 and for the Dirichlet boundary when g > 1 but not for both in the same time. 

The operators involved in the expansions (|5.9| ) and ( |5.17| ) are not the only boundary 
operators presented in the theory. Each kind of boundary allows its special boundary 
operators. 

Consider first the Dirichlet boundary. Since all points have the same x-coordinate, it is 
kinematically impossible to introduce an order operator exp(ipx). However, we can define 
a disorder operator x^m] with magnetic charge m representing a discontinuity Ax = rmx 
at some point of the Dirichlet boundary. Geometrically this operator is represented by a 
source of m domain lines starting at the same point at the boundary. The expectation 
value of such an operator can be calculated by decomposing the world sheet along the m 
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lines (fig. H): 



/>oo 

Ml... I d£m+i 

(5.20) 



/•oo />oo 

dio / d£i... d£m+i 

'0 

-00 r_ poo 



dzi dzm+i lToa[v{-zi, n)]lin[v{-Z2, n)]...lin[v{-Zm+i, /u)] 



M Jm {z + Zi){zi + Z2)...{Zm + Zm+l){Zm+l z) 

Since v{z, n) ~ the whole integral scales as 2;*^"^+^^^/^"^. On the other hand, the 

amplitude with a marked point on the Dirichlet boundary scales as dv{z, i~i)/dz ~ 
Comparing the two powers we find the dimension of the desorder operator on the boundary 

5[m] = mg/2 (5.21) 

With the Neumann boundary the things stay in the opposite way. The disorder 
operators do not make sense because there is already a discontinuity at each point of the 
boundary. However, the order operator V(p) introducing a factor exp[i{p— ^Pq)x] at some 
point ^ of the boundary can be defined perfectly well. Going to the Fourier space and 
distributing the exponential factor among the domain lines crossing the way between the 
point ^ and the Dirichlet boundary, we arrive at the following geometrical description of 
the order operator with electric charge (momentum) p. The expectation value 

i)(p)(z,^) = (V(p)C(^)0(z)) (5.22) 

is equal to the statistical sum for the mixed Dirichlet-Neumann loop amplitude with the 
fugacity of some of the domain lines modified. Namely, all domain lines surrounding the 
point ^ are taken with a factor cos(7rp)/ cos{^TTpo). (We remind that the factor of cos(|7rpo) 
has been absorbed in T) 

The loop amplitude ( |5.22| ) satisfies the following integral equation 

V(^p){z,^)= -J r . p7 (5.23) 

Jm cosi^npo) ■k{z + z') [v{z,h)\^ 

which can derived in the same way as eq. (|3.8|). Eq. ( |5.23| ) can be considered as the 
dispersion integral for an analytic function with a cut M < z < 00. Therefore along the 
cut we have 

n - ( M cos(7rp) lui[w{z)]v^p){-z,n) 

Im^;(p) 2,// = — , z<M 5.24 



Inserting the relation 



lm[v{z, fi)] + ^77^^ = 5.25 
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in (|]2D we find 



-r r-. / XT cosfTrp) Im.\v(z, u)] ^ , , 

Im i)(p)02,;u = -f-^— , z < -M 5.26 



For z large n) ~ z^/^ and therefore 

lm.\v{z,^)] 1 /In /r 

^ sin(-7rfif) = cos(-7rpo), 2^-00 (5.27) 

Therefore at large z the amplitude ( |5.22| ) satisfies the functional equation 

Im[i)(p) {z, n)] = cos(7rp)t)(p) {-z, /j) (5.28) 

whose solution is any power z^ with a = ±{p — ^) + even integer. The leading power at 
^ 00 large can be fixed by the requirement that when the momentum p coincides with 
the background momentum ^po, the amplitude ( |5.22|) coincides with the expectation value 
of the identity operator —d/dn 

d 

t'(ip„)(2, ^) = -^^(^' ^) (5-29) 

Therefore V(p-j {z, n) behaves for z large as 

7)(p)(-;z,;u) ~^|P-^Pol-f (5.30) 

Comparing this with the asymptotics dv/dn ~ ;zS/2-3 we find 

\p-Pq/2\ 

S{p) = 5.31 

9 

For finite z ( |5.22|) is given by the infinite series 

= A)^l^-^^o|-.(H'=)- (5.32) 

k,n 



6. Open string propagator and the spectrum of momenta 

All amplitudes involving closed and open strings are defined by imposing appropriate 
Dirichlet and Neumann boundary conditions on the boundaries of a world sheet with given 
topology. 

Any string amplitude can be decomposed into elementary pieces (propagators and 
vertices) following the logic of refs. [|15[ and . 



In this paper we concentrate ourselfs on the calculation of the open string propagator. 
It will be obtained following the same steps as in the case of the closed string propagator 
T5[] . Contrary to our expectations, the case of open strings turned out to be technically 
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much more difficult tlian tlie case of closed strings. We have found the spectrum of the 
propagator but we were not able to obtain the explicit form of the eigenstates. 

Before considering the open string, let us repeat the major steps of the calculation of 
the closed string propagator, using the SOS model. One has to calculate a string-string 
amplitude with the world sheet configuration of the loops as shown in fig. |^, with non- 
contractable domain walls going around the cylindric surface. In this way we take into 
account the possibility for the closed string to propagate in the x-space. 

Let X and x' be the coordinates of the two (Dirichlet) boundaries of the cylinder. If 
we pass to the momentum space the factor e^'^pi^-^ ) can be written as a product of factors 
g±i7rp associated with the domain walls wrapping the cylinder. Taking into account the 
two different orientations, each such domain wall acquires a factor 2cos(7rp). Further, the 
amplitude of each elementary cylinder between two subsequent noncontractable domain 
walls is [0 



M(£+£') 

G,{lA') = Vl (6.1) 

This amplitude describes the deformation of the closed string from the " in" state of length 
I to the "out" state with a length without a change of the x-space position. The whole 
propagator G{j)\ £, £') in the momentum space can be obtained by sewing such elementary 
cylinders: 



G{p;iJ') = V / ... / -^...^[2cos{np)rGo{i,ii)Go{ii,i2)...Go{in,i') (6.2) 

71=0"^^ ^1 

To calculate it we have to diagonalize Go{i,i') in the £-space. This was done in 

r°° 1 

Go{iJ')== dE{i\E) {E\i') (6.3) 

Jo 2cos(7r£/J 



where 



{e\E) = -^TvEsmh{TvE) K,e{M£) ^{my^,e^O (6.4) 

TV 



form a complete set of delta-function normalized eigenstates 

I j{e\E){E'\i) = 2nS{E,E') (6.5) 

It is convenient to introduce the Liouville variable cp = log i; then the integration measure 
becomes uniform: 

df 

j = d(P; l = e't' (6.6) 

The wave functions behave as plane waves in the limit £ {(f) =—>■ — oo) and decay 
rapidly when (p ~ log(l/M). Therefore the 5— function is produced only by the small-£ 
behavior and the normalization coefficient is not affected by the form of the eigenstates for 
i ~ 1/M. This important feature of the half-space quantum mechanics was emphasised 



and well explained in |25]. It can help to calculate directly the spectrum of the kernel Go 
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from its small-£ behaviour. Indeed, let us expand the r.h.s. of (6J-) in Ijl' , assuming that 



I < I' . One finds, writing the series as the result of a contour integration 



Once the irreducible part Go is diagonalized, the r.h.s. of ( p?^ ) can be evaluated 
immediately 

/•oo ^ 

G{p-iJ')= dE{i\E) r^{E\i') (6.8) 

Jq cosh(7r£/J — cos(7rpJ 

The quantum number E plays the role of the momentum of an additional "Liouville" 
dimension. 

The propagator ( |6.8|) is universal in the sense that it does not depend on the back- 
ground momentum; it is the same for any closed string theory with effective dimension 
(central charge of the matter in the language of 2d gravity) less than 1. 

The poles of the propagator define the possible Liouville energies corresponding to 
given momentum p 

iE^ip) = ±p + 2k, A; e Z (6.9) 



If we consider a theory with a background momentum po = \g — 1\ the allowed momenta 
are p = npo, n G 2. For each value p of the momentum together with the lowest energy 
states E — ±p (the creation and annihilation operators of a closed string "tachyon" ) 

there is an infinite discrete set of states which describe infinitesimal deformations of the 
boundary at this point. A boundary of finite length can be expanded as an infinite series 
of such operators. 

Now let us calculate the propagator for the open string. A typical configuration of 
the loops on the world sheet is shown in fig. ^ The generic loop configuration involves 
three kinds of domain walls: closed loops, lines ending at the same boundary and lines 
connecting two different boundaries of the world sheet. This last kind of domain walls 
describes the propagation of the open string in x-space. If we denote by ro{i,i') the 
amplitude of propagation between two consequent domain walls with the lengths £ and £', 
then the full propagator is given by 




(6.10) 



i.e., by the same expression as ( |6.2| ) , with Go replaced by ro/2 cos(7rpo/2) and the measure 
d£/i replaced by d£. The difference between the two measures of integration is due to the 
cyclic symmetry of the closed string which is absent for the open string. Now let us 
calculate To (£,£')• This is the loop amplitude for a disk with a boundary divided into 
four segments having alternatively Dirichlet and Neumann boundary conditions. Such an 
amplitude can be decomposed as a convolution of an amplitude with Dirichlet boundary 
and a number of amplitudes with Dirichlet-Neumann boundaries. The decomposition can 
be performed by cutting the world surface along the most internal open lines. Summing 
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over the numbers m and n of such hues at the two opposite boundaries we find 

k,m=0^^ i=l 
m km 

II dije-^''"'^V{ij)W{i + £' + Y,^^ + Y.^^^ 

j = l i=l j = l 



(6.11) 



oo 



m 



k=0 



dt ^(£+£')w^+^rT>/o^.- ,Mfc+mr2 cos(7rpo/2) . 



f:^e'^'+''^W{t)[V{2Ko-t)y 
^ — ^ / z7^^ 



2ni [1-2 cos(7rpo/2)y(2i^o - t)/T]^ 

where product from 1 to is assumed 1. The contour of integration goes around the cut 
[^Lj^i?] of Wo{t), leaving outside the cut [2Ko — tii,2Ko — ti] of the denominator in the 
integrand. Note that Fq depends only on the sum i + i' of its arguments. 

In the scaling limit we replace the quantities in the integrand by their singular parts 
according to eqs. and ( ^.2tj| ). The regular part in the denominator cancels and we 

obtain the following scaling limit of the irreducible part of the string propagator 

— —To{i,i')= dz e-^^'+''^To{z) (6.12) 

2cos(7rpo/2) Jm 

T. / N 1 lToaw{-z) 1 lm[v{-z,^)] 

^o[z) = - — = — ^ 1, z^oo (6.13 

2cos(7rpo/2j [^(-2jJ 2cos(7rpo/2j v{z, fi) 

(The asymptotic value at z ^ oo follows from the large z asymptotics of the open string 
background amplitude v{z) ~ z^^'^. 

Before proceeding further let us notice that at small lengths ro(^ + i') is identical to 
its analog (|6.1D for the closed string, up to the factor y/IF 

-1 -M{e+e') 

-ro{iJ')= , , „ (6.14) 



2cos(7rpo/2) ' ' ' £ + £' 

The factor ^/JF will appear if we replace the open string integration measure di = e'^dcj) 
with dljl = d(j). Hence one can expect that the spectrum of the open string propagator 
is the same as the one of the closed string, since the spectrum is defined by the small I 
asymptotics of the propagator. However, the eigenfunctions will be certainly different. 

A rigorous proof of this statement can be done by demonstrating that the traces of 
all powers of the two propagators are the same. 

Let us consider the trace of the open string propagator (|6.10| ) and perform the £ 
integrations using the form ( 6.12 ) of Fq. The result is 



/"OO 

Trr(p) = / diG{p]i,l 

J a 

= V[2cos(7rp)]'' / ... / dzi...dz 



ro(.i)...ro(..) ^^-^^^ 



(Zi + Z2){Z2 + Zs)...{Zn + Zi] 
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The integral is logarithmically divergent for large z and has to be cut off at z ~ 1/a. 
If we plug the large z expansion 



ro(2) = l + (const.)2"^ + ... (6.16) 

in ( |6.15| ) , the constant term reduces to exactly the same convolutive integral as in the case 
of the closed string propagator. This term diverges as log(l/a). All subdominant terms 
will produce jinite corrections. Therefore 

Trr(p) = lim log(-) / — - + finite terms (6-17) 

a^o a Jq cosh(7ri?) — cos(7rp) 

One can repeat the same argument for any integer power of the pro pag ator. This means 
that the open string propagator can be written in a form similar to (|6.8[) 



POO ^ 

r{p;ij')= / dE{i\E)A,^ {E\i')A,>. (6.18) 

Jq cosh{TvE) — cos{iTp) 

All dependence on the cosmological constant A and the mass /x at the ends of the string is 
absorbed in the eigenstates of the propagator 

{E\i)A,f, = {V,EC{i^)0{i))A (6.19) 

Of course, the next terms of the asymptotics also contain some universal information 
about the open string theory, and they will be important for the calculation of the string 
interactions. Note that the on-mass-shell (microscopic) states E = ±ip are the order 
parameter amplitudes ( |5.22| ). 

Once the propagator is known, the amplitudes involving string interactions can be 
calculated by decomposing the world sheet into irreducible pieces (vertices and propaga- 
tors) in the same way as it has been done in the case of the closed string. The dependence 
on the external momenta is only through the propagators. For the three-string amplitude 
this is illustrated by fig. ^ The vertices rn(^i, ^i, ^2, ^2, ••, ^n+2, ^n+2) represent amplitudes 

for a disk with n + 2 pairs of Neuman and Dirichlet boundaries with lengths £k and ik, 
/c = l,2,...,n-|-2. Itis easy to see that these vertices will depend only on the total lengths 

£ = ii + ... + in+2 and i = ii + ... + in+2 of the Dirichlet and Neumann boundaries. 
Therefore the Laplace image of (we denote it by the same letter) will depend only on 
two variables z and fu,. It is easy to establish the following integral equation which follows 
from geometrical decomposition of the world sheet shown in fig. ^ 



Therefore 



Imtf) ( z) 

ImTniz, a) = — ... along the cut z < M (6.21) 

[v{-z,^)]''+^ 



In particular, the tadpole vertex {n = —1) is the basic open string amplitude 

F_i(z,^) =0(2,^) (6.22) 
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and the integral representation ( |6.20|) becomes a closed equation which is the continuum 
limit of {^). 

The spectrum of excitations of the open string is fixed by the set of allowed target 
space momenta p . Since the background momentum po/2 of the open string is twice less 
than that of the closed string, its spectrum will be twice denser 

P = o^Po, n G 2 (6.23) 



7. Conclusion 

We have demonstrated in this paper an approach to the open string theory which 
allow us not only to calculate the scaling dimensions of the boundary operators, but also 
to obtain, in principle, any given amplitude for the open strings in the dimensions less than 
1, with arbitrary in- and out- momenta. The propagator of two open strings presented 
here is the simplest possible example. Technically, this question is not simple, since one 
has to know the eigenfunctions of this propagator, for which we know only the functional 
equation. 

Another interesting possibility is to try to solve the whole field theory for the open 
strings, which seems now to be a difficult but not a hopeless task. This might shed some 
light on the string picture in the multicoulour QCD according to an observation made in 
1?^ . This might be an object of further study. 

Finally, let us note that the spectrum of the open string excitations is not related to 
the value of the parameter fj, which can be considered as the mass of the "quarks" at the 
ends of the open string. 
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Appendix A. 

To give a more precise meaning to the derivation of the eq. (|3.2| )for the open string 
amplitude let us use its definition by means of a three-positioned graphs as shown in the 
fig. |TI|. A world sheet of the string amplitude Kri,™, looks like a (p^ planar Feynman diagram 
with rh legs occupied by the ends of open lines (thick lines) and m nonoccupied legs ( thin 
lines) at the boundary, rh and m are the lengths of the Neumann and Dirichlet boundaries, 
correspondingly. To every link on the graph occupied by a loop or open line one subscribes 
a weight 1/(2X0)- 

If we pick up one leg at an edge of the Neumann boundary, the corresponding open line 
decomposes a graph of the amplitude into two similar ones, only with the different lengths 
of the Neumann and Dirichlet boundaries, and this geometrically obvious decomposition 
allows us to write the loop equation in the form: 

m 00 f' I V 

VrH,m = J2Y1 Wk-2,,WrH-k,p+J^^^{2Ko)-P-'^-^ (A.l) 
k=2p,q=0 
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If we introduce the generating function 



-771 — 1 



m,m=0 



m=2 



m=0 



p,q=0 



plql 



(A.2) 



this equation transforms in the following way 

^ oo CX3 ^ r J r 



(A.3) 



where Vfc(t) is the amplitude with k legs on the Neumann boundary {k is even, of course) 
as a function of the spectral parameter t of the Dirichlet boundary. In particular, Vo(t) = 

W{t) is the amplitude with pure Dirichlet boundary. Performing all the sums and the 
integration over a in ( A.3 ), we arrive to the following equation 



TViT, t) = Wit) + I ^ J-y(T, r)V{T, 2Ko - r) 

J ZlTl t — T 

which is identical (up to a normalization of V(T, t)) to our main equation ( p.5|) . 



(A.4) 
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Figure Captions 

Fig. 1. The geometry of the world sheet with alternative Dirichlet-Neumann boundary 
conditions. The thick lines represent the Neumann boundaries 

Fig. 2. Critical points of the embedding of the world sheet, a) Creation (annihilation) 
of a closed string state, b) Splitting (joining) of closed strings, c) Creation 
(annihilation) of open string state, d) Splitting (joining) of open strings. 

Fig. 3. A configuration of domain walls for a world sheet with Dirichlet-Neumann bound- 
ary conditions 

Fig. 4. The geometry of the loop equation for the Neumann-Dirichlet disk. The fat line 
represents the Neumann boundary. 

Fig. 5. The geometry of the classical field equation for the open string tadpole 

Fig. 6. The geometrical description of a disorder operator and the corresponding decom- 
position of the world surface 

Fig. 7. Loop configuration for the closed string propagator 

Fig. 8. A typical configuration of domain walls for the open string propagator 

Fig. 9. Decomposition of the world sheet for the three string amplitude 

Fig. 10. The decomposition of the world sheet producing the integral representation for 
the vertex Fi 

Fig. 11. Discretized world sheet of the open string amplitude in the loop gas representation 
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